Effect of mesoscopic inhomogeneities on local tunnelling density of states 
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We carry out a theoretical analysis of the momentum dependence of the Fourier-transformed local 
density of states (LDOS) in the superconducting cuprates within a model considering the interference 
of quasiparticles scattering on quenched impurities. The impurities introduce an external scattering 
potential, which is either nearly local in space or it can acquire a substantial momentum dependence 
due to a possible strong momentum dependence of the electronic screening near a charge modulation 
instability. The key new effect that we introduce is an additional mesoscopic disorder aiming to 
reproduce the inhomogeneities experimentally observed in scanning tunnelling microscopy. The 
crucial effect of this mesoscopic disorder is to give rise to point-like spectroscopic features, to be 
contrasted with the curve-like shape of the spectra previously calculated within the interfering- 
quasiparticle schemes. It is also found that stripe-like charge modulations play a relevant role to 
correctly reproduce all the spectral features of the experiments. 

PACS numbers: 74.25.Jb, 74.20.-z, 74.50.+r, 74.72.-h 



I. INTRODUCTION 



In the last few years scanning tunneling microscopy 
(STM) measurements have become a most valuable tool 
in investigating the physical properties of superconduct- 
ing cuprates. Although, like angle-resolved photoemis- 
sion spectroscopy (ARPES), this technique is mostly 
sensitive to the surfaces, the layered structure of the 
cuprates suggests that it may be representative of the 
bulk properties. In particular, in the recent years Fourier- 
transformed (FT) STM spectra showed a rich momen- 
tum structure in the local density of states (LDOS) and 
raised a strong debate on what elementary excitations 
produce such structure a 1 : 2 : 3 ^ 5 : 6 . To be specific the dot- 
like patterns observed in FT-STM spectra have been at- 
tributed either to local charge-spin order with pinned 
static collective excitations^ or to interference effects be- 
tween impurity-scattered quasiparticles (QP) 1 ' 2 : 5 . Sev- 
eral theoretical analyses have tried so far to consider 
these mechanisms and to relate them to the experimental 
observations. Grossly speaking, it is found that pinned 
collective textures may account for the rather punctual 
(although obviously broadened) character of the LDOS 
patterns in momentum space 7 ' 8 ' 9 ' 10,11 . However, the pre- 
dicted patterns show too weak an energy dependence, 
which contrasts with the substantial dispersion of most 
of the experimentally detected spots. On the other hand, 
the theoretical analyses based on QP interference gener- 
ically produce dispersive LDOS patterns, but at a given 
energy the high intensity regions form extended curves 
in k spaceiSiiiiiii^, which hardly resemble the experi- 
mental spot-like intensity patterns. Moreover the STM 
measured dispersion curves do not properly match the 
ARPES-determined QP dispersions and tend to "flatten" 
at wave- vectors typical of charge/spin order. Finally, the 
weight of QP peaks in ARPES is weak and strongly de- 
pends on temperature. Therefore the QP interference 



effects should disappear upon approaching T c 1517 . This 
seems to be the case for many of the structures, but ex- 
perimentally it is also found that some of the k-space 
features persists above T<&. Therefore neither of the two 
pictures is fully satisfactory and one still needs to rec- 
oncile these analyses with experiments. In this paper we 
precisely aim to perform a systematic analysis at low tem- 
perature of the QP picture to determine whether or not 
physically sensible mechanisms (atomic form factors, dis- 
order, multiple scattering, mesoscopic inhomogeneities, 
charge-ordering instabilities) can turn the curve-like QP 
spectra into more spot-like patterns. We consider the 
structure factors due to the short-distance structure of 
the Wannier orbitals, the second-order impurity scatter- 
ing processes, the effect of magnetic impurities and, most 
importantly, the effect of mesoscopic inhomogeneities. In 
this way we succeed in smoothing most of the curve-like 
LDOS features into broad peaky structures. However, a 
detailed comparison of the experimental figures with our 
calculated spectra shows that some features cannot be 
properly reproduced. A better agreement is instead ob- 
tained when additional effects from charge-density mod- 
ulations are considered. Therefore, our analysis shows 
that in the real cuprate systems there must be a coex- 
istence between dispersive QPs (producing dispersive in- 
terference patterns broadened by inhomogeneity effects) 
and incipient static local charge order responsible for the 
enhancement of non-dispersive peaks in specific regions 
of the k space. 



The paper is organized as follows. Section II contains 
the model and the description of the approach. In Sec- 
tion III the mesoscopic inhomogeneities are introduced 
in the model and their effects are described. In Section 
IV we present our results, while our concluding remarks 
are contained in Section V. 
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II. THE MODEL AND THE TECHNIQUE 
A. The model 

In STM experiments, the LDOS is measured with 
atomic resolution on the points r of a field of view I x I 
with t ~ 640 A (several hundreds of lattice unit cells). 
Once the LDOS N(r, u>) is obtained, its Fourier transform 
is a function of the momenta (q x ,Qy) = {n x ,n y )2Tr/£, 
yielding the wave- vector power spectrum 



P(q,u;) = \N(q,co)\ 2 /e 2 



(1) 



Notice that, since the r positions of the STM scans are 
denser than the atomic positions, the q momenta are not 
restricted to the first Brillouin zone. In the following we 
work on an infinite lattice defined on atomic positions R 
(in units of the lattice spacing a) and use a field of view 
of size L x L, where L — t/a. Therefore the momenta 
are restricted to the first Brillouin zone. This restriction 
is relaxed in Appendix A by including the atomic form 
factors encoding the orbital subatomic structure. 

To consider the elastic scattering of the quasiparticles 
on quenched impurities, we introduce a (weak) external 
potential e(R) = e (R- — R«) due to the local poten- 

tial of Ni impurities randomly located on sites R^ of the 
L x L — N two-dimensional system. Although most of 
the expressions below stay valid for a general form of the 
impurity potential, to be more specific we consider the 
form 



e(R) = J2 V 8(R - Ri) 



which in momentum space reads 



e(q)=^oE eiqR ' 



(2) 



(3) 



Aiming to perturbatively calculate the corrections in- 
duced by the impurity potential on the LDOS, we need 
to calculate the electron Green's function in the super- 
conducting state. Thus it is convenient to introduce 
fcrmionic Nambu spinors to write the Green's functions 
in matrix form 

0(k,w) ^(k,w) 

where the normal and anomalous Green functions are 
given by 



(4) 



Also the scalar (i.e., non magnetic) impurity potential 
can be put in matrix form 



e s (q) = e(q) 



allowing one to define the quantity 
d 2 k 



A_(q,u/) = J 



(2n) 2 



(5) 



G(k, u/)£(k + q, u) - F{k, w).F(k + q, u) 
d 2 k (oj + e k )(w + £k+ q ) - A k A k+q 



(6) 



(2n) 2 (^-E^-E 2 ^) 
as the (1,1) element of the matrix product 

g(k,u>) ^(k, w ) 

1 0\f £(k + q,w) .F(k + q,w) 
-1 J \T(k + q,uj) £(k + q,-w 

The integral in Eq. (JSJ emphasizes the fact that we are 
working on an infinite lattice. 

Since the LDOS is obtained from the imaginary part 
of the (1,1) element of the Green function matrix, at 
first order in the impurity potential, the correction to 
the LDOS can be written as 

d 2 q 



iV (1) (R,w) 
d 2 q 



(2n) 2 



=> R e(q)A_(q, W ) 



(2tt)2 

{(£ K (q) cos(q • R) + e 9 (q) sin(q • R)) 9f(A_(q,«)) 
+ (e^q) sm(q • R) + e S (q) cos(q • R)) K(A_(q,cj))} 
Moreover, since 



(7) 



e^( q )= [e(q) + e(q)*] = -[e(q)+e(-q)] 



2 v ^ v ^ J 2 
1 „ 1 



(8) 



e w (q) = - [e(q) - e(q)*] = - [e(q) - e(-q)] (9) 

then e 5R (q) is the symmetric part of e(q) with respect 
to q — > q while e 9 (q) is the antisymmetric one. Both 
9?(A_(q)) and 5R(A_(q)) are symmetric with respect to 
the same transformation so the only term that survives 
under the sum over q is the first one in Q, then 

(e s (q)cos(q ■ R) + e 9 (q) sin(q ■ R)) 3(A_(q,w)) 

If there is inversion symmetry (for example one symmet- 
ric impurity at R = (0,0)) so that e(R) = e(— R) then 
£ J (q) — an( i e ( ( l) — e5R (q) an d one simply has 

d 2 q 



jV (1) (R, 



{2-k) 2 

For Ni delta-like impurities one has 



e(q)cos(q.R)3(A_(q,^)) (11) 



Ni 



c R (q) = VoX) c08 ( < l- R <) 

i=i 

e 9 (q) = K.2sin(q-R i ) 



(12) 
(13) 



i=l 
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where are random positions. 

It is worth noticing that in this scheme the LDOS 
N(R,u>) is calculated for a fixed configuration of disor- 
der (i.e. of impurities) and no average over these config- 
urations is taken. Furthermore we neglect the contribu- 
tions from impurities that are outside the field of view. 
This results in Fourier-transformed scattering potentials 
e(q), which could be sizably (and randomly) momentum- 
dependent. Only owing to the rather large size of the 
L x L field of view and to the self-averaging character 
of this disordered system, the e(q) functions turn out to 
be sufficiently smooth to preserve the momentum struc- 
ture of the 3(A_(q, oj)) encoding the interference effects 
of the scattered QP's. In particular for an infinite sys- 
tem one would obtain \e(q)\ 2 /N w (\e(q)\ 2 )/N = mV^ 
(where the angular brackets denote disorder average )i»V 
One could also calculate the LDOS power spectrum 
by Fourier transforming the LDOS correlation function 
(N(R,uj)N(0)) . Again, this procedure would only allow 
to extract informations on the momentum structure of 
A(q) if carried out on sufficiently large grids such that 
the resulting e(q) is smooth. 

We consider a system with a bare tight-binding band 
ffe = — 2i(cos k x + cos k y ) — 4i' cos k x cos k y — fi (t and t' 
are the nearest-neighbor and the next-nearest-neighbor 
hopping parameters on a square lattice and fiQ is the 
chemical potential) and a <i-wave superconducting gap 
Afe = Ao(cosfc;c — cosky)/2. In this way the QP disper- 
sion is i?k = \/ e k + ^k- To compare our results with 
the experiments of Refs. [U l'Jol we take the following pa- 
rameter values, which, for t w 150 meV are suitable for 
a Bi2Ca2SrCu20g+x sample around optimal doping 16 : 
t' = -0M, A = 0.25t = 37.5 meV. For the sake of 
definitencss we take the doping x = 0.15, for which we 
calculate the chemical potential to be = — l.Ot. 

Fig. 1 displays the momentum dependence of the 
LDOS for this system and well reproduces typical re- 
sults of Ref. Capriotti et ali^. To make the comparison 
with figures reported in experimental papers, throughout 
this paper we will mark as darker the regions with larger 
spectral intensity. Moreover, we orientate the momen- 
tum axes in such a way that the Cu-0 directions (usu- 
ally taken as the x = (1,0) and y = (0, 1) directions in 
theoretical papers) here are taken along the diagonals of 
the figures. Notice also that, with the chosen orientation, 
the second Brillouin zone is visible in our figures. 

The spectrum in Fig. 1 has been discussed in terms of 
the constant-energy curves of the quasiparticle dispersion 
relation E^. Fixing the E^ identifies "banana" shaped 
curves in momentum space with large density of states 
at the extremities^*^. 

These spectra are quite different in various respects 
from the experimental ones. It is apparent that the 
high- intensity regions are extended lines, which do not 
reproduce the spot-like shape of the experimental high- 
intensity regions. One can show that large contributions 
to the spectra are obtained when both denominators in 
Eq. [S] are small. Extended one-dimensional features are 




( - ih , ) ( -7t , -Jt) ( , - 2 71 ) 



FIG. 1: Momentum-dependent LDOS for an homogeneous 
Bi2Ca2SrCu2 08+x superconducting system with x = 0.15. 
The momenta are in units of the inverse lattice spacing. To 
facilitate the comparison with experiments in this and subse- 
quent figures we rotate the Brillouin zone so that the Cu-O 
directions of the square CuC>2 planes are along the diago- 
nals of the figure. The tight-binding parameters are (t = 
150 meV) t' = -0.3t; no = -l.Ot; A = 0.25t = 37.5 meV, 
u — — 0.08i = —12meV. The concentration of impurities is 
1%. 



obtained when a translation of q of one banana makes 
it tangent to another banana. Changing q keeping the 
two bananas tangent defines a one-dimensional feature of 
high intensity. Crossing of two of these features produces 
a high intensity spot. However those spots are quite dif- 
ferent from the experimental ones since they are clearly 
associated with the one-dimensional crossings. 

Another important difference is that there is no inten- 
sity around zero momentum. This lack of spectral weight 
in the q = (0, 0) region contrasts with the presence of 
rather intense broad peaks appearing in the experimen- 
tal data. 

Finally, the theoretical treatment considers a mesh co- 
incident with the atomic positions and correctly repro- 
duces a spectrum, which is periodic in momentum space: 
N(q, u>) = A^(q+G, oS) with G a reciprocal lattice vector. 
However, this feature is not present in the experimental 
spectra, where the peaks in the second Brillouin zone 
are suppressed with respect to their first Brillouin zone 
partners. The presence of these form factors is easily ac- 
counted for by the local space structure of the Wannier 
orbitalsii and this effect is described in Appendix A. As 
shown in Fig. 2, the introduction of Wannier orbitals 
does not significantly improve the unrealistic appearance 
of the spectrum. This is most evident at low momenta, 
where the non-point-like nature of the orbitals is obvi- 
ously immaterial. 

In the next Sections, we will elaborate on the above 
expressions to consider the inhomogeneous distribution of 
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( 0,2%) I 

I 



( -2^, 0) (-%,-%) (0,-2Tl) 



FIG. 2: Momentum-dependent LDOS spectra for a system 
with the same parameters as in Fig. 1 and Gaussian Wan- 
nier orbitals with width a — 0.0625 in units of square lattice 
spacing. The impurity concentration is 1% and the spectrum 
is taken at uj — — 0.08i = — 12meV. 



doping and the consequent inhomogeneous distribution 
of the chemical potential and of the superconducting gap. 



III. EFFECT OF MESOSCOPIC 
INHOMOGENEITIES 

STM experiments show large fluctuations in the gap 
amplitude over large length scalesiS*^. It is quite impor- 
tant to recognize that the size of these regions is of several 
unit cells, £ ~ 13 — 15^4, and is apparently unrelated with 
the average distance between impurities. 

Competition among different phases in strongly 
correlated systems can give rise to mesoscopic 
inhomogeneities2ii22i2i. One can expect that this 
effect is enhanced by the inhomogeneous distribution of 
the doping. Here we will not discuss the microscopic 
mechanisms that can give rise to this effect but consider 
it as granted and we analyze the consequence on the 
LDOS. 

Since the mesoscopic inhomogeneities involve several 
unit cells and differences can be substantial from one 
region to the other we choose an approach, which is 
analogous to the usual semiclassical treatment of elec- 
trons in presence of slowly varying perturbations, where 
the electron distribution depends both on space and mo- 
menta. In the same spirit we allow the Green functions 
to depend parametrically on the real-space region of the 
sample via the space dependence of both the chemical 
potential /i = /i(R) and the superconducting (SC) gap 
A = A(R). This gives rise to a local density of states 
which depends explicitly on R due to the conventional 
impurity scattering, and implicitly trough the paramet- 



ric dependence of the gap and the chemical potential. 

tf(q,w) = ^e-^ R iV(A(R), A i(R),R,^ (14) 

R 

Here the sum is restricted to the L x L field of view and 
the first order expression of iV(A(R), /i(R), q, oj) is given 
by Eq. [7| but with the Green functions computed with 
the local value of /i(R) and A = A(R). For example Q 
is given by Eq.(0J, but with 

€k = — 2t(cosk x + cos ky ) — 4t'cosk x cosk y — ^i(R) 

(15) 

Afe = A(R)(cos k x — cos k y )/2 

These expressions introduce a parametric dependence of 
Q on the position R via the local values of the chemi- 
cal potential and of the maximum value of the d-wave 
SC gap Ao (however, to keep the notations simple, in 
the following we often do not explicitly indicate this R 
dependence). This approach allows for substantial differ- 
ences in the gap and local chemical potential from one 
region to the other and hence go beyond conventional 
perturbative formulations. 

To realize a specific inhomogeneous distribution of 
fi and A we consider space fluctuations of the doping 
around a given average value x. Specifically we generate 
an inhomogeneous map characterized by doping fluctu- 
ations of 30% (i.e. x = x(R) locally ranges from 0.1 
to 0.2) of typical size £ ~ 3 — 4 lattice units. We take 
this range of fluctuations as an estimate deduced from 
the (larger) relative fluctuations of thegap observed in 
the Bi2Ca2SrCu208+ x sample of Ref. [l8| and under the 
assumption that the gap and the doping are linearly re- 
lated. For simplicity the smooth map was replaced by a 
mesa like function by determining contour levels of the 
smooth map and assigning to all points between two suc- 
cessive contour levels a constant doping equal to the av- 
erage value of the two limiting contours. In this way 
the doping interval was (arbitrary) coarse-grained in five 
slices and for each of the five possible doping values the 
corresponding values of the chemical potential and of the 
SC gap where calculated. In particular, starting from the 
given tight-binding structure, the chemical potential was 
determined according to the local doping x(R), while for 
simplicity the maximum value of the d-wave gap was de- 
termined by a linear rescaling with doping. The resulting 
coarse-grained map is shown in Fig. 3 for a 200 x 200 field 
of view. A direct comparison shows a close resemblance 
between our space inhomogeneity map and the similar 
experimental figures of Refs. llSllliJ . 

A. Zeroth order 

Even in the absence of impurity scattering, the inho- 
mogeneities on the chemical potential and on the SC gap 
affect the LDOS spectra. In particular if £ is the typi- 
cal size of the domains in which fi and A can be taken as 
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nearly constant, it is quite natural to expect that a broad- 
ening of the spectra is obtained around the q — (0, 0) 
wave- vector, over a k-space range of the order of £ . 
Specifically, 

2V (0) (q,w) = ^e- iqR 7V (0) (A(R),^(R),cj) (16) 

R 

7V(°)(A(R),^(R),w) is the LDOS without impurities for 
a (homogeneous) system with SC gap A(R) and chemical 
potential £t(R): 



AT(°)(A(R), M (R)) 



d 2 q 
(2tt) 2 



g(A(R) lM (R),q,o;) 



where the Green function is taken in the absence of im- 
purity scattering and is given by the same expression of 
Eq.Q, but with space-dependent parameters /x(R) and 
A(R). 

The numerical evaluation of JV' ' (q, lj) [Eq. I|lt)|) ] for 
the specific distribution of inhomogeneities shown in Fig. 
3 (upper panel), indeed produces a rather strong peak 
around zero momentum (Fig. 3, lower panel). This ef- 
fect can be simply understood as follows. The inhomo- 
geneities in A(R) and /i(R) reflect in the DOS, which 
deviates from its average value Nq. For illustrative pur- 
pose, we consider an approximate linear dependence of 
the DOS from the energy in the nodal approximation, so 
that, for small deviations of A and /x we find 



7V(A(R), M (R))=7VoM 



1 



A(R)-A , M (R) 



Mo 



/'o 



(17) 

Here Ao = (A(R)) and /io = (M^-)) arc the average val- 
ues of A and /i respectively. We also assume a Gaussian 
distribution 



(18) 



where A is suitably chosen to give the expected size of 
the fluctuations of about thirty per cent. Here X(R) = 
<5A(R), fyt(R), with SA(R) = A(R) - A and Sn(R) = 
/i(R) — /io. The LDOS correlation function then takes 
the form 
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FIG. 3: Upper panel: Real-space map of a 200 x 200 field of 
view with randomly distributed regions of tipical size £ ~ 3.5 
(in units of lattice spacing) and five different values of dop- 
ing and corresponding values of gap and chemical potential 
(see text). Lower panel: Momentum-dependent LDOS with- 
out impurity scattering as it arises from Eq. 1161 with the 
mesoscopic distribution of Fig. 3 (upper panel). 



(JV(A(R),/i(R))JV(A(0), M (0))> = 



(19) 



1 



(<5A(R)c5A(0)) (tf/i(R)5/z(0)) 



A 2 



l4 



Taking the Fourier transform, one obtains the LDOS 
power spectrum 



P(q,w) = JV (w) 



A 



A 



en 2 



i 



1 



(20) 



Therefore the [i and A distributions generate the peak 
around q = (0, 0) in the figures reported in the next 
subsections. As expected, we find that the width of this 
peak scales as £ -1 showing that this simple effect explains 
the central bump, which is not reproduced in the usual 
"homogeneous" approaches. 



B. First order 

Besides the above simple zero-th order effect, the local 
inhomogeneities affect the first-order corrections induced 
by the scattering of the QP on the random impurities. 
Therefore, in the presence of the mesoscopic inhomo- 
geneities, we recalculate the first-order corrections to the 
LDOS 

NW(H,U>) = -I^ e ^ R 3 [ / ^e^ R e (q')A_(q> 

(21) 

where A_(q, ui) has the same form as in Eq. JSJ), but 
with independent normal and anomalous Green func- 
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tions according to their parametric dependence on A(R) 
and /i(R). The sum over R is to be intended over the 
points of the L x L grid of the field of view. We also 
consider the possibility of QP's scattering on magnetic 
impurities, for which the matrix representation is given 
by 



£m(q) = e m (q) 



1 
1 



In the case of magnetic impurity scattering, the expres- 
sion of A for the LDOS of a spin-up electron reads 



A+(q,w)= J 



d 2 k 



G R (k, uj)G R (k + q, w) + ^(k, w)-F fl (k + q, oj) 
d 2 k (uj + e k ){u + £k+ q ) + A k A k+q 
(2n) 2 {^-El){^-El +ti ) 



(22) 



Since an overall minus sign is obtained for spin-down elec- 
trons, no contribution is obtained at first order for the 
total LDOS in the case of magnetic impurities. 



C. Second order 

At second order in the impurity concentration, if we 
consider only scattering processes that occur on the same 
site, we have 

N^U^lu) = -- Ve- !ClR x (23) 

71 R 

^ {/ (2^ elPR Pl (P ' ^ ^ (P ' ^ + S2(P ' W) A+ (P ' " )] } 



where 



Exfow) = ( e 2 (p)+e 2 m (p))^ 

k 

E 2 (p,w) = (e 2 (p) + 4(p))£ 



uj 2 - El 



(24) 



(25) 



and A± are given by Eq.© and Eq. respectively. 

Various observations are in order here. First of all 
it is found that at second order, a finite contribution 
to the LDOS in momentum space is obtained also in 
the case of magnetic scattering. Notice also that, con- 
trary to the first-order case, the real parts of A_ and 
A+ both contribute to the LDOS. Finally we notice from 
Eqs. (|23|) - (|25[) that second-order QP interference pro- 
cesses may contribute to the LDOS at zero momentum 
N(q = (0,0)). This contribution adds to the zeroth- 
order peak arising from the mesoscopic inhomogencitics 
described in Sec. III.B and may be present even in the 
absence of such inhomogeneities. This can be seen by 
taking /i(R) = an d A(R) = A . 

Although at q = (0,0), the first-order contribution to 
A_(0) is small [see Fig. 4 (upper panel) and Appendix B], 



A+(0) does not vanish and contributes to A^ 2 ^(0) (both 
in the absence and in the presence of magnetic scatter- 
ing). Therefore second-order scattering processes could 
contribute to the intense (rather broad) peak experimen- 
tally obtained for N(q,uj) at q = (0,0). However, we 
checked from the relative intensity of the peaks at finite 
momenta that the scattering processes are weak and the 
second-order processes are not strong enough to explain 
the rather large intensity of the peak at q = (0, 0). This 
strongly indicates that the zero-order contribution from 
the mesoscopic inhomogeneities is the main source of the 
large intensity of the zero-momentum peaks in the exper- 
iments. 



IV. MICROSCOPIC SPATIAL CHARGE 
MODULATIONS 

The observation (see below) that calculations of the 
STM spectra in terms of interfering scattered quasipar- 
ticles do not account for the observed weight of some 
specific spectral features motivated a further enrichment 
of our treatment. Specifically we considered a physi- 
cally different effect arising from microscopic charge mod- 
ulations. Since long time it has been suggested that 
the superconducting cuprates might be close to an in- 
stability leading to the spatial ordering 24 likely in the 
form of fluctuating onedimcnsional charge textures (the 
so-called stripes), which can even acquire a slow criti- 
cal dynamics 2 ^ around optimal doping and above. The 
scale of these textures is of a few lattice spacing thus we 
call them "microscopic" as opposed to the mesoscopic 
inhomogeneities consider above. The proximity to the 
microscopic instability can naturally reflect itself in a 
large enhancement of the charge susceptibility at spe- 
cific wave-vectors Q c ^ (obviously corresponding to the 
charge order) and to related structures in the momentum- 
dependent dielectric function. In particular, it is natural 
that the external potential introduced by static impuri- 
ties is substantially screened by the nearly unstable elec- 
tron liquid and acquires some structure in its momentum 
dependence. One can see that close to the charge insta- 
bility, the static charge susceptibility acquires a nearly 
polar form2& 



x(q) = 



xo(q) 



xo(q) 



l-V(q)Xo(q) |q-Q c fc| 2 +C/f 



(26) 



where £ c /j is the correlation length for charge fluctuations. 
Once this charge susceptibility is introduced in the di- 
electric constant, one finds that the impurity potential is 
screened as 



Vo(q) = V 1 + 



1 



ch) 



(27) 



1 and 



For simplicity, we assume here VoXo(q = Q 
we take = 0.25 inverse lattice spacing. This screening 
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gives rise to a strongly momentum-dependent impurity 
potential. Quite obviously, once this screened impurity 
potential is inserted in the expressions for iV(q, us), it will 
filter the momentum dependence of the LDOS emphasiz- 
ing the intensity at q « Q c h- 



V. RESULTS 

We here describe the effects in the calculated 
momentum-dependent LDOS arising from the progres- 
sive introduction of the mesoscopic inhomogeneities and 
of the spatial charge modulations according to the scheme 
described in Section III and IV. To compare with the re- 
sults reported in Figs. 1 and 2 for an homogeneous sys- 
tem, we keep considering the tight-binding model with 
the same parameters. Again the value of no corresponds 
to a doping x — 0.15. Besides the previously considered 
effect of the local Wannier structure of the orbitals, we 
introduce a random inhomogencous (mesoscopic) doping 
distribution according to the scheme of Subsection III. A. 
Fig. 4 represents a momentum-space LDOS including 
up to second-order scattering. As it can be clearly seen, 
the space inhomogeneity is quite effective in modifying 
the STM spectra (a) by introducing substantial spectral 
weight around zero momenta, and (b) by broadening the 
curve-like spectral features. In particular, these latters 
acquire a "fuzzier" appearance, which emphasizes the re- 
gions of stronger intensity thereby producing more spot- 
like features and rendering the calculated spectra closer 
to the experimental data. We also notice that this effect 
arises from the disorder in gap and chemical potential 
and is therefore physically quite different from (and its 
effects quite more pronounced than) the impurity disor- 
der considered in the appendix of Ref. flH 

Despite this substantial improvement, a closer com- 
parison with the data of Refs. and y (cf. Fig. 5 
(lower panel) shows that the calculations of Fig. 4 fail 
in reproducing some of the rather intense features exper- 
imentally detected. Specifically, strong spectral features 
are observed along the Cu-O-Cu directions (i.e. along 
the diagonals of Figs. 4 and 5, see the experimental in- 
tensities schematically reported in Fig. 5, lower panel) 
at wave-vectors corresponding to a four-unit-cell mod- 
ulation |q| = 0.25 x (27r/ao) (ao is the lattice unit in 
the supposedly square CuOi planes) 4,5-2 -, which are to 
weak in Fig. 4 (lower panel). Therefore, following the 
scheme of Section IV, we consider the effects of screen- 
ing on the impurity potential. According to the sugges- 
tion that the cuprates are close to a charge instability, 
this screening is taken to be strongly momentum depen- 
dent, as a result of the strong charge susceptibility of 
the electron liquid at momenta q rj Q c h- The resulting 
calculated spectrum is reported in Fig. 5 (upper panel), 
where the too weak peaks in the (±1, 0), (0, ±1) direc- 
tions (along the diagonals of the figures) are now strongly 
enhanced. The comparison with experimental intensi- 
ties is more direct in Fig. 5 (lower panel), where we use 



( 0, 2n) ) 



f 



(-2k, 0) (-71, -Jt) (0,-2%) 



(0,271) ■ 



( -71 » 7t) 




(-2~k,0) (-7C,-7l) (0,-27t) 



FIG. 4: Upper panel: First-order-only contribution to the 
momentum-dependent LDOS spectra for a system with the 
same parameters as in Fig. 1 and 2 (Wannier orbitals are 
also considered). The chemical potential and the SC gap are 
randomly distributed according to the real-space map of Fig. 
3 (upper panel). The impurity concentration is 1% and the 
spectrum is taken at uj — — 0.08i = — 12meV. Lower panel: 
Momentum-dependent LDOS spectra with the impurity scat- 
tering calculated including the zeroth, first and second order 
with e = 0.8i. All other parameters are the same as in the 
upper panel, 

dots to depict the regions where experimental spectra dis- 
play the strongest features. It is clear that the effect of 
momentum-dependent screening on the impurity poten- 
tial is substantial and important to correctly reproduce 
all the experimental structures. 

VI. CONCLUSIONS 

The results displayed in the previous sections allow 
to draw some general conclusions concerning the mecha- 
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( 0, 2n) ). 



( -it, t) 




(-21,0) (-7t,-7l) (0,-2Jl) 



( 0, 2Jt) L 



( -7C/ ft) 




(-2*1, 0) (-7t,-7C) (0,-27l) 



FIG. 5: Upper panel: Momentum-dependent LDOS spectrum 
including the effects of mesoscopic inhomogeneities and the 
impurity scattering up to second order with e — 0.8t. We also 
considered a momentum-dependent screened impurity poten- 
tial according to Eqs. 12711 and . 1271 with f = 0.25 in units 
of inverse lattice spacing. The impurity concentration is 1% 
and the spectrum is taken at ui = — 0.08t = — 12meV. Lower 
panel: same as in the upper panel, but with the addition of 
dots marking intense features in the experiments of Ref. 

nisms leading to the formation of (point-like) structures 
in the LDOS as obtained in Fourier-transformed STM 
spectra. 

First of all, disregarding the zero-order inhomogenity 
effects acting at low momenta, we considered the inten- 
sity of the finite-momentum peaks. From the comparison 
of the first and second-order contributions with experi- 
mental spectra, we can draw the conclusion that the QP 
scattering due to impurities is rather weak. Therefore 
our perturbative approach is appropriate and the first- 
order processes already account well for the part of the 
spectra, which can be attributed to the QP interference. 
When this mechanism at first order fails in reproducing 



the all relevant features of the spectra, the second-order 
processes do not substantially improve the calculations. 
The same holds true as far as the inclusion of magnetic 
impurities is concerned. We also took in consideration 
the extended character of the Wannier orbitals, which 
partially modifies the spectra at large wave-vectors ac- 
counting for the weakening of the spectral features at 
large momenta. Nevertheless the calculated spectra pre- 
serve their unrealistic features. 

A substantial improvement in the theoretical calcu- 
lation is represented by the inclusion of space inhomo- 
geneities both in the SC gap and in the chemical poten- 
tial as arising from doping inhomogeneities. As far as 
the impurity scattering is concerned, the main effect of 
inhomogeneities is to blur the curve-shaped spectral fea- 
tures promoting a more point-like appearance of the high 
intensity regions of the spectra. This is quite easily inter- 
preted in terms of the (by now) standard arguments re- 
lated to the relevance of some specific wave- vectors join- 
ing the parts of the electronic (d-wave gapped)spectra 
at fixed energy bias (see, e.g., Fig. 1 of Ref. .5). The 
contours of constant-energy in the quasiparticle spectra 
have the well-known "banana" shape and in the absence 
of disorder the momentum space regions with large in- 
tensity are extended curves. However, by introducing 
local gap fluctuations, the length of the "bananas" fluc- 
tuates, while the transversal, i.e. the width, fluctua- 
tions of the "bananas" are instead less pronounced due 
to the very elongated shape of the constant-energy con- 
tours. In this case one sees that the gap fluctuations 
produce already a marked blurring of the LDOS peaks. 
The introduction of the local chemical potential fluc- 
tuations, instead, tends to produce shifts in the Fermi 
surface and in the quasiparticle spectra, which are most 
pronounced along the Fermi velocities vp . This locally 
induces transversal shifts of the "bananas" , further blur- 
ring the LDOS spectral features. As a result, the unreal- 
istic curve-like regions of high LDOS intensity in momen- 
tum space are transformed in spot-like features, which 
have a closer appearance to the experimental spectra. 
This is the first main outcome of our work. Despite this 
success, however, there are some specific spectral fea- 
tures at some specific momenta, |q| = 0.25 x (27r/ao) 
(close to the so-called qi in the "octet" model of Ref. 
y) which we did not reproduce with the proper inten- 
sity within the QP interpherence mechanism. It was al- 
ready observed in Ref. 0, that theoretical calculations 
based on interference effects between impurity-scattered 
quasiparticles underestimate the above spectral features. 
Moreover, recent experiments^ in the pseudogap region 
of Bi2Ca2SrCu208+ x above have shown the persistence 
of LDOS peaks at some momenta, which can hardly be 
attributed to interphering QP's and could rather be at- 
tributed to some spatial order in the charge and/or in the 
spin channels. Our findings strengthen this observation 
showing that the several improvements at large (Wannier 
functions) and small (zero-order effects) wave- vectors, as 
well as second-order calculations do not change this con- 
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elusion. Therefore, aiming to reproduce all the features 
of the STM spectra, we included a momentum-dependent 
screening due to a (supposedly strong) charge suscepti- 
bility at some specific wave- vectors. As described in Fig. 
5 above, the close resemblance of the calculated spectra 
to the experiments suggests that a strong tendency of the 
cuprates to order spatially at some specific wave- vectors 
might well be present in the real systems and coexist 
with d-wave quasiparticles. This important indications 
represents the second main outcome of our work. 

After this work was completed, we became aware of 
more recent STM experiments 2 ^ where two different spa- 
tially separated regions were identified in underdoped 
Bi 2 Ca2SrCu 2 08+x samples. In (more metallic) regions 
with large coherence peaks (which are the large major- 
ity around optimal doping) the dispersive feature related 
to the qi vector of the "octet" model has a smooth in- 
tensity variation. On the other hand, in (less metallic) 
regions where the coherence peaks are absent, the disper- 
sive spectral feature related to qi acquires an additional 
sizable intensity, when qi is close to the incommensurate 
charge ordering momenta q w 0.25 x (27r/ao)(±l, 0) and 
q w 0.25 x (27r/a )(0, ±1). This suggest a substantial 
charge ordering in these regions. From this point of view 
our Fig. 5 represents a "superposition" of the spectra 
from these spatially separated regions. 
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APPENDIX A: LOCAL DENSITY OF STATES 
WITH WANNIER FUNCTIONS 

The first striking difference with experiments is that 
A_(q, to) is periodic in the reciprocal space, whereas the 
experimental spectra are not. One first natural step to- 
ward a more realistic description of the STM spectra can 
be carried out by considering the real space structure of 
the Wannier orbitals. In this way the short-distance de- 
scription around each lattice point is improved, leading to 
a refinement of the calculated STM spectra at large wave- 
vectors. Defining 0R,(r) a Wannier function with center 
in R, a vector of the lattice, the electron Green functions 
can be expanded on this basis so that the LDOS at the 
first order in the impurity concentration can be written 



in the most general way 

(r,w) = 3 J dr G(r,ro)e(ro)G(r ,r, W ) (Al) 
= 3 E E E {G(k 1; c))G(k 2 ,a;))x 

ki,k 2 Ri,R 2 R'j ,R 2 

J dro #R.i ( r o) e(ro) <pn 2 (r )| 

where GG = QQ in the normal phase and GG = QQ — TT 
in the superconducting phase If one assumes that the 
impurity potential e(r) is only defined on the lattice sites 
R, then 

f dr^,(r)e(r)0 R2 (r) ~ J RiiRa £ e"^ e(q) (A2) 
J q 

and one obtains the following expression for the LDOS 
N(t,u) = E G(k,^)G(k + q,c)e lk - R (A3) 

k,q R-! ,R.;> 
e - 4 (k +q ).R 2£(q) ^ Ri(r) ^ (r) 

If one further assumes that the Wannier orbitals are 
rather localized around each lattice site, the following 
approximation becomes applicable 

^RiW^W^^R^I^RiWI 3 (A4) 

and taking advantage of the property of the Wannier 
functions 

&i(r)=&i+R'(r + R'), (A5) 
the Fourier transformation of N(r, u), 

N(q,w) = J dre^ r N(r,uj) (A6) 
can be written as 

N(q, u>) =F(q)J2 e~ <qR 9 E ^ A - (<*> ^ A7 ) 
R q 

Here the form factor 

F(q) = J dre- tcl r |</» (r)| 2 ~ e -^ +q y^ (A8) 

has been made explicit by assuming a Gaussian shape 
of the (rather localized) Wannier orbital, with a propor- 
tional to the spacial variance of the Wannier functions. 

It is quite simple to check that, if Wannier functions 
are delta-like, a — * and F(q) — > 1, and the usual ex- 
pressions of Ref. are recovered. 
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APPENDIX B: LOCAL DENSITY OF STATES IN 
THE POINT (0,0) 

We numerically find [Fig. 4 (upper panel)] that 
A_(q = 0, to) is quite small. This numerical finding can 
acquire an analytic support, from a calculation at low to, 
where a nodal approximation is justified. In this case we 
can show that A_(q = 0, w) vanishes. Let us consider 



A±(q,w) = 



dk (uj + ek)(uj + 6k+ q ) T A k A k+q 



{ ^-El){^-El + J 



At first order only A_ enters the calculation for non- 
magnetic impurities, while only A + is to be used for 
magnetic impurities. In this latter case the summation 
over spins leads to a cancelation and no contribution is 
obtained at first order. On the other hand, at second or- 
der, both A_ and A+ are needed both for magnetic and 
non-magnetic impurities [cf. Eqs. (|23|l - i|25j . In order to 
calculate the local density of states at q = (0, 0) both at 
first and second order we consider 



A±(0,w) 



dk (w + t k ) 2 TAj 



(2tt)2 ( W a-2£) a 
If we change the variables of integration near a node 



(Bl) 



dk 



2 71 



dB 



dpp 



(2ir) J 2~KVpv\ 



so that 



£k = P cosy 
Ak = p sin#, 



the Eq. I|B1(I becomes 

A±(q,w) = 4 / — x 
Jo ( 27r ) 

r dp 2 ^+p 2 ((cosg) 2 T(sing) 2 ( ^ ) 
J Attvfva (w 2 — p 2 ) 2 



For A_(0) we have 



A-(0) 



27r de 

M W) 
1 



dp 2 



p 2 cos 20 



ttvfva 



dp 2 



u 2 

{uj 2 - p 2 ) 2 



p 2 f 



(B3) 



For A+(0) we have instead 



A+(0) = 



ttvfva 



dp 2 



(lu 2 - p 2 ) 2 



The local density of states in the long wavelength limit 
in q = (0, 0) and in the presence of non magnetic and 
magnetic impurity is given by 



TTVpVA 
1 



dz- 



dz 



(lu 2 — zy 
uj 2 + z 



= 



1 



(u> 2 — z) 2 VfVA 



(B4) 



(B5) 



■KVpVA 

with a double pole in uj 2 + i0 + . This means that the 
spatial average of n(R) at first order is zero for a non 
magnetic impurity, while is finite for a magnetic impurity 
at fixed spin. In this latter case it could provide a direct 
measure at low energies of the product between Fermi 
and gap velocities. However, as already stated at the 
end of Sect. III.C, summing over spins one again obtains 
a vanishing contribution. On the other hand, at second 
order both A_ and A + contribute to the scattering in the 
presence of magnetic and/or non-magnetic scattering. 
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